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ABSTRACT 

We demonstrate that observations of glitches in the Vela pulsar can be used to 
investigate the strength of the crust-core coupling in a neutron star, and suggest that 
recovery from the glitch is dominated by torque exerted by the re-coupling of superfluid 
components of the core that were decoupled from the crust during the glitch. Assuming 
that the recoupling is mediated by mutual friction between the superfluid neutrons and 
the charged components of the core, we use the observed magnitudes and timescales 
of the shortest timescale components of the recoveries from two recent glitches in 
the Vela pulsar to infer the fraction of the core that is coupled to the crust during 
the glitch, and hence spun up by the glitch event. Then, within the framework of a 
two-fluid hydrodynamic model of glitches, we analyze whether crustal neutrons alone 
are sufficient to drive the glitch activity observed in the Vela pulsar. The analysis is 
conducted using two sets of neutron star equations of state (EOSs), both of which 
span crust and core consistently and cover a conservative range of the slope of the 
symmetry energy at saturation density 30 < L < 120 MeV, a range which encompasses 
all current experimental constraints. One set produces maximum masses ~2.QMq, the 
second ~2.QMq, allowing us to probe the effect of the high-density stiffness of the 
EoS as well. We also include the effects of entrainment of crustal neutrons by the 
superfluid lattice. We find that for medium to stiff EOSs, observations imply > 70% 
of the moment of inertia of the core is coupled to the crust during the glitch, though for 
softer EOSs L « 30MeV as little as 5% could be coupled. No EOS is able to reproduce 
the observed glitch activity with crust neutrons alone, but extending the region where 
superfluid vortices are strongly pinned into the core by densities as little as 0.016fm“^ 
above the crust-core transition density restores agreement with the observed glitch 
activity. 

Key words: dense matter - equation of state - starsmeutron - pulsars: general - 
pulsars: individuahVela 


1 INTRODUCTION 


Glitches - sudden increases in the rotational frequency of 
pulsars - offer an insight into the internal dynamics of neu¬ 
tron stars and potential constraints on the properties of 
dense matter. The quasi-periodic giant glitches exhibited by 
the Vela pulsar ( Melatos et al.|2008 1 offer some of the most 
stringent tests of our pulsar glitch models, and our mod¬ 


els of dense matter (Link et al. 19991. A leading class of 


glitch model supposes a component of the stellar interior 
spends most of the time decoupled from the solid crust and 
magnetosphere whose secular spin-down we track observa- 
tionally (see [Haskell fc Melatos (20151 for a review of glitch 
models). This component, usually taken to be part of the 
neutron superfluid, occasionally recouples to the crust and 


transfers some of its angular momentum, spinning the crust 
up ( [Anderson fc ltoh[1975||Alpar|1977| . The superfluid com¬ 
ponent can be decoupled by the pinning of vortices to nuclei 


in the inner crust (Alpar 

1977 Pines et al.||1980 Anderson 

jet al.[[1982| [Alpar et al.| 

1984P, or by pinning to magnetic 

flux-tubes in the core (Ruderman et al. 1998 Link 20031. 


As a rotational lag builds between the decoupled compo¬ 
nent and the rest of the star, so the hydrodynamic Magnus 
force builds until it is able to overcome the pinning force and 
cause the vortices to move outwards collectively, transferring 
angular momentum to the crust. Due to the complexity of 
physical glitch scenario, only recently have significant steps 
been made in a consistent hydrodynamical simulation of the 
glitch process ([Andersson et ar]|2006[ [Sidery et al.|[^10[ 
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van Eysden fc Melatos|[2010| [Haskell et al.|[^12[ [Haskell ^ 

Antonopoulou|2014 1. 


In the most widely examined variant of the model, it 
is the crustal superfluid that is pinned. However, recent cal¬ 
culations which show that a large fraction of the superfiuid 
neutrons are entrained by the crust via Bragg scattering off 
the crustal lattice, has led to a vigorous debate on the ef- 
hcacy of this mechanism in the light of the observed Vela 


glitch activity ( 

Charnel 

|MT3 

2012 

Andersson et al. 

2012 

Piekarewicz et al.[[2014 

[Steiner et al.[[2015[ [Hooker et al.[ 

2015 

). From this, one can conclude that if the whole star 


spins up during the glitch, it appears that the existence of 
entrainment in the crust renders the crust-driven glitch sce¬ 
nario only marginally viable at best. If, however, a portion of 
the core is decoupled from the crust during the glitch, there 
might be sufficient angular moment stored in the crust su¬ 


perfiuid to drive the observed giant glitches (Hooker et al. 

[Ms] ). 

It is expected quite generally that the core superfiuid 
will be coupled to the normal (non-superfluid) fluid by the 
mutual friction force. The main mechanism which is thought 
to be acting in the core is scattering of electrons off super- 
fluid vortex cores, magnetized by entrainment of protons. In 
this case Mutual Friction will couple the crust and the core 
superfiuid in less than a minute for a typical glitching pulsar 
with a rotation period of around 10 ms ([Alpar et al.||1984[ 


Alpar & Sauls 1988 Andersson et al. 20061. In the crust. 


however, protons are not superfiuid and electron scattering 
is ineffective. In this case the Mutual Friction coupling is 
likely to be much weaker and due mainly to interactions with 
sound waves in the lattice ( Jones|19'OT I. Vortices close to the 
crust-core boundary, that are mostly immersed in the crust, 
would thus be weakly coupled, with the coupling becoming 
stronger at higher densities where vortices are mostly im¬ 
mersed in the core and electron scattering is the dominant 
Mutual Friction mechanism ( Haskell et al.|2012| ). Compari¬ 
son with the observational upper limit on the timescale for 
Vela glitches <40s ( [Dodson et al.|[2002| ) suggests that it is 
likely a significant part of the core will not be coupled to the 
crust during the glitch. This uncoupled component should 
then recouple in the minutes after the glitch, a process that 
should manifest itself observationally in the recovery of the 


pulsar spin frequency immediately post-glitch (Haskell et al. 


2012 Haskell & Antonopoulou 20141. Indeed, in the 2000 


and 2004 Vela glitch, exponential fits to the recovery of 
the spin frequency provide tentative evidence for a short 
timescale (~ Imin) component ( Dodson et al.|2002 20071. 

Recently, much progress has been made simulating the 
hydrodynamic evolution of the superfiuid vortices includ¬ 
ing the coupling of crust and core via mutual friction and 
microscopic pinning forces. Such hydrodynamic models can 
explain qualitatively the inter-glitch timescales. Vela glitch 
sizes, and post-glitch rotational evolution, and have the po¬ 
tential to constrain the neutron star equation-of-state (EOS) 


(|van Eysden &; Melatos 

20101 |Haskell et aL|2012||Pizzochero| 

2011 Seveso et al.|2012 

1, despite remaining uncertainties in 


aspects of the glitch model such as the unpinning trigger 
mechanism and details of how the vortices subsequently un¬ 


pin (Glampedakis & Andersson 2009[ Warszawski & Melatos 

2008 Melntos_&_WaxszawsM 2009[ [Warszawski fc Melatos 


2013 Warszawski et al. | |2012|lLink||2014[ ), as well as not 

including the effect of crustal entrainment. 


In this paper, we use the qualitative picture that 
emerges from hydrodynamic models together with the short 
timescale component of the 2000 and 2004 Vela glitches to 
infer the moment of inertia of that portion of the core super- 
fluid that is coupled to the crust during the glitch via mutual 
friction. We then calculate the maximum moment of iner¬ 
tia of the crustal superfiuid neutrons in the crustal regions 
which drive the glitch according to hydrodynamic models, 
including the effects of crustal entrainment, and hence in¬ 
fer the glitch activity for the Vela pulsar which we confront 
with the observed value. The neutron star models used are 
generated using two families of equations-of-state; one gen¬ 
erated using the Skyrme non-relativistic nuclear model and 
the second using the relativistic mean field (RMF) model. 
Both families predict the same EOS up to saturation densi¬ 
ties, but differ at high densities; the RMF model is stiff and 
able to accommodate maximum neutron star masses above 
2.5Mq, while the Skyrme model used is softer and gener¬ 
ally gives maximum masses around 2.QMq. We generate the 
families of EOSs by varying the slope of the symmetry en¬ 
ergy at saturation density, a quantity that has been shown to 
correlate strongly with radius and crust-core transition den¬ 
sity. We vary it over a conservative range 30-120MeV which 


encompasses constraints inferred from experiment (Tsang 


et al. 2012 Lattimer & Lim 20131. When L is varied, the low- 


density pure neutron matter is adjusted to maintain a good 
fit to the results of ab initio pure neutron matter (PNM) 
calculations (Gezerlis & Garlson 2010 Hebeler & Schwenk 


2010 Gandolfi et al.||2012 l. In doing so, we are able to ex¬ 

plore the predictions of our glitch model over a wide range 
of EOS parameter space. 

In section 2 we explain how we use the observed recovery 
to set constraints on the fraction of the core coupled to the 
crust at the time of glitch. In section 3, we apply this to our 
glitch model using our consistent sets of EOSs. We present 
our results In section 4, and discuss them in section 5 as well 
as giving our conclusions. 


2 CRUST-CORE COUPLING VIA MUTUAL 
FRICTION 

The Vela pulsar spins with an angular frequency of H « 
70 rad/s, and exhibits giant glitches (AH)g « 1.5 x 10“"* 
rad/s. The time taken for a glitch to occur (the glitch rise 
timescale) is constrained by observation to Tg < 40s. The 
secular spin-down rate due to el ectromagnetic torque on th e 


Vela pulsar is H 


10 - 


(Dodson et al. 


2002 


20071. 


These are taken to be the values for the solid crust of the 
Vela pulsar, to which the magnetic field is anchored. The 
normal component of the star is coupled to the crust on 
short timescales, but the superfiuid neutrons are only weakly 
coupled, and may not react fast enough to remain coupled 
during a glitch. Let us take as an illustrative lower limit only 
1% of the total moment of inertia of the star to be coupled 
to the crust during the glitch. Then, by angular momentum 
conservation, one would expect the immediate post-glitch 
spin-down rate to be tl ~ 10“® s“^. 

The key assumption in our model is that the core super- 
fluid recouples post-glitch via mutual friction, in which su¬ 
perfluid vortices become magnetized by entraining protons, 
allowing electrons to scatter off them and hence couple to 
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Figure 1. The post glitch spin down for the total model (with 
mutual friction and electromagnetic torques included) and for a 
model with only mututal friction included. The spin frequency is 
taken to be that of the Vela pulsar, $1 = 70 rad/s. We see that 
for the first hour or so mutual friction dominates the spin down 
and the two curves are indistinguishable. After approximately an 
hour the electromagnetic torque drives the spin down and the 
curve for the total model begins to show faster spin down. 


the charged components of the core and, through them, the 
crustal lattice. We assume that the protons and electrons 
in the core are rigidly rotating and coupled to the crust on 
short timescales by the magnetic field ( Haskell et al.|[^012 
van Eysden|[2M4| . 


Let us now examine how effectively mntual friction con- 
pies the crnst and core on short timescales. Mutual friction 


couples the two components on a timescale | Andersson et al. 


[2006] ) 


"rmf ' 


1 

206 


( 1 ) 


where B is the mutual friction coefficient (Andersson et al. 


2006) which expresses the strength of the mutual friction 


force. Setting this equal to the glitch rise timescale Tg ~ 
Imin gives us a value of 6 ~ 10“"^. All regions of the core 
with B > 10~'* will be coupled to the crust during the glitch, 
while regions of the core with B < 10“^ will begin recoupling 
to the crust after the glitch, resulting in a torque on the 
crust. A region with B « 10“®, for example, will recouple to 
the crust on a timescale of several minutes after the glitch, 
and the resulting torque will give a contribution to the spin 
down rate, due to mutual friction, of 


^ 2nBAn « 10"'^S“^ 


( 2 ) 


This is considerably stronger than the electromagnetic con¬ 
tribution, even assuming a small amount of the moment of 
inertia of the core is coupled during the glitch, and we can 
can infer that core recoupling via mutual friction will give 
the dominant contribution to the spin-down rate for up to 
several hours after the glitch. 

These estimates are indeed confirmed by running the 
code of|Haskell et al.|(|2012|) for a Vela giant glitch. In Fig 
0 we compare the spin down rate after the glitch obtained 
for a standard simulation and the rate obtained by excluding 
the external electromagnetic torque. We can see that the two 
rates agree in the hrst few hours in which mutual friction 



2000 

2004 

AFg(Hz) 

2.45435x10-® 

2.2865x10-® 

AUtr (Hz) 

2x10-® 

5.4x10-® 

T (min) 

1.2±0.2 

1±0.2 


Table 1. Constant step in frequency Fg and fast decaying com¬ 
ponent of the step Ttr, with associated timescale t, for the Vela 
2000 and 2004 glitches. 


dominates, and only after do they differ as the electromag¬ 
netic contribution takes over. Thus immediately after the 
glitch, the crust component of the star (including the core 
protons) is spinning faster than a portion of the core neutron 
superfluid. 

If this is the case the observation of short timescale spin 
recovery after a glitch represents a direct probe of the mutual 
friction mechanism, and therefore the vortex dynamics in 
the outer core. We shall now demonstrate how one can infer 
the moment of inertia of that part of the star coupled to the 
crust during the glitch, 7c. It is this quantity which enters 
into the parameter G = 7sf/7c (where 7sf is the moment of 
inertia of the superfiuid component driving the glitch) to be 
compared with the pulsar’s glitch activity. 


2.1 Extracting the coupled moment of inertia 
from Vela glitch observations 

There are two published glitches that have an observed 
short-term recovery component, the Vela 2000 and 2004 
glitches ( Dodson et al. |2002[ |2007| ) . Unfortunately, as shown 
in table the amplitudes of these components, ATtr vary 
signihcantly between the two by more than a factor 1000, 
partially due to the fact that in the Vela 2004 the short term 
component was barely above the noise level. The timescales 
r are, however, similar, and consistent with each other 
within the errors. 

Let In be the moment of inertia of the superfluid re¬ 
gion that is in the process of recouping on the timescale Tmf 
(~ Imin from table [^, providing the observed short-term 
contribution to the spin-down of the star, while 7c is the 
moment of inertia of the crust and all components tightly 
coupled to it during said timescale (which includes the glitch 


event itself). These are related through (Andersson et al. 

[2006l l 


2116 (l-f 7„/7c) • 

We shall refer to these two components as n and c. 

Before the glitch, the two components will have a small 
equilibrium lag (AS7)eq = « 10“® for the Vela pulsar. 

The glitch acts to spin up the component c by an amount 
(Af2)g 3> (Af2)eq, and therefore immediately after the glitch 
the lag between the components will be (Afl)in = (All)g. 
After the timescale over which the component n recouples, 
component c has decreased its frequency by an amount 
(Af2)tr and component n has increased its frequency by 
(Af2)g — (Afl)tr; thus by angular momentum conservation 
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/n ^ (An)tr 

7 c ~ (An)g - (An)tr 

Taking (Af2)g and (Af2)tr from observations, we obtain 
^ « 8 X 10”^ for the Vela 2000 glitch 

J-C 

^ « 2.37 for the Vela 2004 glitch 

-^C 


(4) 

(5) 

( 6 ) 


Using this result and the observed timescales, equa- 
tion[^gives us values for the mutual friction parameter above 
which correspond to regions of the core coupled at the time 
of glitch 


75 = 1x10 ^±3x10 ® for the Vela 2000 glitch (7) 

B = 4x lO”'^ ± 1 X 10"® for the Vela 2004 glitch (8) 

These numbers are still (roughly) consistent given the 
uncertainties that are clearly underestimated in this proce¬ 
dure. In the following we shall simply assume a range: 

3 X 10"® < Bohs < 1.3 X 10“'‘. (9) 

from which we can compute the moment of inertia of the 
core 7c coupled via mutual friction to the crust at the time 
of glitch, and hence spun up by the glitch. This procedure, 
although limited by current observational uncertainties, al¬ 
lows us to derive this important quantity from observations 
for the hrst time. 

We now outline how we use this result to extract 7c 
from neutron star models. 


3 MODELING THE GLITCH 

Given a particular neutron star equation of state (EOS), we 
obtain our background neutron star model by solving the 
Oppenheimer-Volkoff (OV) equations, thus obtaining the 
density and composition profile of the star as a function 
of radial coordinate r. 

The moment of inertia of a star of radius R in the limit 
of small angular frequency Q (Hartle & Thorne 19681 is given 
by 


-^tot - 


3 Jo ^ - 2GM{r)/r ’ 


( 10 ) 


where e{r) is energy density of matter in the star, P{r) is 
the pressure and M{r) is the mass contained in radius r. 
i'{r) is a radially-dependent metric function given by 




2GM\ H {M{x) + ‘inx^P{x)) 
~ir )~^ Jr *2 (1 _ 2GM{x)/x) 

( 11 ) 


and uj is the frame dragging angular velocity 



Figure 2. Neutron star cross-section in plane of rotation axis (Q) 
depicting the geometry of the strong pinning region in inner crust 
(shaded area SP) and the region of the core coupled to the crust at 
the time of the glitch via mutual friction (shaded region MF). OC 
and IC label the outer and inner crust respectively. The angular 
locations of the boundary of the strong pinning region at the outer 
crust, and the outer and inner boundaries of the coupled region 
of the core, are given by Pouter, 0i and O 2 respectively. Except for 
the softest equations of state at saturation density (L ps 30 MeV), 
92 = 0 . 


for r ^ R. 

We take the standard form for the mutual friction coef¬ 
ficient in the core, due to electron scattering off magnetized 
vortex cores ( Andersson et al.|2006 1: 

2 / \ 1/2 

.. . . / m_ \ ' 

B = 4 X 10” 


/^W/6 / 

V 0 . 05 / V 


V lO^'^g/cm® 


(14) 


where Xp is the proton fraction, m* the microscopic effective 
proton mass and p the total density in the core. These quan¬ 
tities are obtained consistently from the EOS, as a function 
of radial coordinate r once the OV equations are solved. 


As detailed in Haskell et al. (2012|, we average B over the 


z-direction to find the average mutual hction strength expe¬ 
rienced by a vortex at cylindrical radius f = rsin0. 


B{f) = 


f 

J 7 


2Z3 


rdd 


S(r) _ 

^ 1 — 2M (r)/r 

^ =rd6 


(15) 


We can now define a region of the star bounded by 
cylindrical radii fi and f 2 such that B{ri < f < r 2 ) > 
7?obs. This corresponds to the region in which the superfluid 
neutrons are strongly coupled to the crust on the glitch rise 
time. 

Dehning 


i_ A 

r3 



dcu(r) \ 
dr ) 


+ 4 


dj(r) 

dr 


uj(r) = 0, 


where 


( 12 ) 


(13) 


^ 4 -,.(r) Cu(r) {e„{r) + P„{r)) 

3 ® H y/l-2GM{r)/r' 


(16) 


where Sn{r) is the energy density of the superfluid neutrons, 
Pnir) is the pressure of the superfiuid neutrons, then the 
moment of inertia contained in the regions of the star with 


j{r) = e = \/l — 2GM{r)/re 


































Crust-core coupling during glitches 5 


cylindrical radius greater than fi and f 2 respectively is given 
by 


c'kI’I pRcc 

Ii ,2 = / / ClsinOdOdr (17) 

-' 01,2 JR(e) 

and therefore the moment of inertia of the region in which 
the superfluid neutrons are strongly coupled to the crust on 
the glitch rise time is given by 


Inc ~ I 2 — Ii- (18) 

Here, Rcc is the radius of the crust-core boundary, and 61,62 
are the angular locations where the outer and inner bound¬ 
aries of the coupled region where it meets the crust-core 
boundary respectively (see Fig. 2). The total moment of in¬ 
ertia of that part of the star strongly coupled to the crust 
at glitch rise time is obtained by adding the contribution of 
the crustal lattice itself and the protons in the core 


Ic — Inc Ic 


. + In- 


(19) 


Hydrodynamic simulations of vortex evolution suggest 
that only the crustal superfluid neutrons within the strong 
pinning region of the crust, defined as the region within 
which vortices are totally immersed in the inner crust, con¬ 
tribute to the glitch itself. How many contribute depends 


on details of the pinning force throughout the crust (Haskell 
et al.|2012| ), but here we take the entirety of the strong pin¬ 
ning region as an upper limit. The moment of inertia of the 
strong pinning region of inner crust superfluid neutrons is 


3.1 Nuclear matter parameters and crust and 
core equations of state 

The glitch model requires several microscopic quantities as 
input. These include the total pressure and energy den¬ 
sity P(nb), £(nb), neutron pressures and energy densities 
Th(nb), en(n-b), effective proton mass in the core, proton 
fraction and mass density mp(nb), a;p(nb) and p(nb), as a 
function of baryon number density nb, and the crust-core 
transition density and pressure rice and Pec- These quanti¬ 
ties are derived from an underlying model of nuclear matter. 

Experimental information about nuclear matter is pre¬ 
dominantly extracted from nuclear systems at densities 
around nuclear saturation density no = 0.16 fm~® and at 
proton fractions close to one half. As a consequence, nuclear 
matter models are generally characterized by their behav¬ 
ior in that region of parameter space. Denoting the energy 
per particle of nuclear matter around saturation density by 
E{nh, S) where <5 = 1 — 2a;p is the isospin asymmetry param¬ 
eter; 5 = 0 corresponds to symmetric nuclear matter (SNM), 
and 5 = 1 to pure neutron matter (PNM). The nuclear sym¬ 
metry energy S{n) is defined as the quadratic coefficient in 
the expansion of E{n\^,S) about 5 = 0: 

Nuclear matter models can be characterized by their 
behavior around nuclear saturation density no = 0.16 fm“®, 
the density region from which much of our experimental in¬ 
formation is extracted. We can denote the energy per parti¬ 
cle of nuclear matter around saturation density by E{nh, 5), 
where nb is the baryon density and 5 = 1 — 2a: the isospin 
asymmetry, where x is the proton fraction, x = 0.5,5 = 0 
corresponds to symmetric nuclear matter (SNM), and x = 
0,5 = 1 to pure neutron matter (PNM). By expanding 
E(n,x) about 5 = 0 we can define the symmetry energy 
S{n), 


r(®p) _ 

■*csf 


nTT/2 j- 
I 0outer I R 


R(0o 

R(8) 


r) 


r Tdr sinM^ 


( 20 ) 


where R( 0 ) is the distance from the core of the star to the 
inner boundary of the strong pinning region at an angle ff to 
the rotation axis, 7?(^inner) ^ Rinner aud 7?(^outer) ^ 77outer 

(see Fig. 2). 

Entrainment of superfluid neutrons by the crust’s lat¬ 
tice reduces the mobility of the neutrons with respect to that 
lattice. It can be shown that this effect is encoded by intro¬ 
ducing an effective mesoscopic neutron mass m n ^ amel| 


2005 Charnel & Carter 2006 Charnel 20121; larger values 


correspond to stronger coupling between the neutron su¬ 
perfluid and the crust, and a reduction in the fraction of 
superfluid neutrons able to store angular momentum for the 
glitch event. One can include this effect by modifying the 
integrand Eq. f^in the inner crust: 


2 _ 

r X ^ r X = 


2 *(r) 


r X 


( 21 ) 


where m*(r) is the effective mass at radius r in the crust. 
We obtain m’^{r) from the results of Charnel ( Chamel|2012 1 
by interpolating between the values calculated at specific 
densities to find the effective mass at arbitrary locations in 
the inner crust. 


E{n, 5) = Eo{nh) -I S{nb)S^ -I- . 


( 22 ) 


The symmetry energy is the energy cost of increasing the 
isospin asymmetry of matter and is a function of baryon 
density. Eurthermore, we can expand the symmetry energy 
about saturation density using the density parameter x = 

'n-b-’T-O - 
3no 


S{TLh) — J Lx hEsynnX “f •■■5 


(23) 


where J, L and Agym are the symmetry energy, its slope and 
its curvature at saturation density. Over the past decade, 
vigorous effort has been devoted to experimentally con¬ 
straining J and particularly L ( Li et al)||2008 Tsang et al. 
20121. Currently the congruence of experimental results 
( Hebeler et al.|2013| ) points to the range 30 < L < 60 MeV, 
but stiffen (higher values of L) are not conclusively ruled 
out ( Fattoyev &; Piekarewicz|2013| ); we therefore examine a 
conservative range of 30 < L < 120 MeV in this paper. 

Crust and core EOSs and the transition density are cal¬ 
culated consistently using two models of the nuclear many- 
body interaction. We use the lUFSU parameterization of the 


relativistic mean-field model (Fattoyev et al. 20101 and a 


parameterization of the non-relativistic Skyrme model, Skl- 
UFSU, which shares the same saturation density symmetric 
nuclear matter (SNM) properties as lUFSU, used in previ¬ 
ous work ( [Fattoyev et al.||2012[ |2013| [Hooker et al.|[^015p . 
Both models have isovector nuclear matter parameters ob- 
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Figure 3. (Color online) Crust thickness ARc (a) and Total radius R (b) versus mass, for the L = 30 (blue), 60 (pink) and 90MeV (blue) 
members of the non-relativistic SklUFSU EOS which is softer at high densities (solid lines) and the relativistic lUFSU EOS (dashed 
lines) which is stiffer at high densities. 



r (km) 



r (km) 


Figure 4. (Color online) Mutual friction profiles as a function of radial coordinate for a 1.4 Mq (a) and 1.8Mq (b) star for the L=30,60,90 
and 120 MeV members of the SklUFSU EOS family. We infer from observations the value above which the superfiuid neutrons are tightly 
coupled to the crust; the shaded region depicts the upper and lower bounds of that value. 


tained from a fit to state-of-the-art PNM calculations (Gez- 


erlis fc Carlson||2010| [Hebeler fc Schwenk||20l0{ [Gandolfi 
et al.|20'T2|, which makes them particularly suitable for de¬ 


scribing the low-density neutron fluid in the inner crust, and 
both describe the bulk properties of doubly magic nuclei 
well ( Fattoyev et al.|[2012 l. The presence in both models of 
two purely isovector model parameters allows the density 
slope of the symmetry energy L to be systematically ad¬ 
justed, while retaining the fit to PNM at low densities. The 
fact that they are isovector means that such adjustments 
leave SNM properties unchanged (Chen et al. 20091. This 
adjustment under the PNM constraint leads to the relation 
J = 0.167L + 23.33 MeV. 


Note that you generally cannot separate out the indi¬ 
vidual proton and neutron pressures in these nuclear matter 
models - the pressure contains terms that are not separa¬ 
ble into neutron and proton components. However, in both 
Skyrme and RMF models, the enthalpy density h = e + P 
is separable, because the non-separable terms in the pres¬ 


sure and energy density cancel out. It is the enthalpy that 
appears in the moment of inertia integrals. 


Both nuclear matter models give closely identical EOSs 
up to saturation density for a given value of L. The SklUFSU 
model is softer than lUFSU at high densities and but gives 
maximum masses of M « 2Mq for all values of L, matching 
the observational lower limit ( [Demorest et al.]|2010|| |Anto-| 
niadis et al.|2013b The lUFSU model we have chosen to be 


maximally stiff at high densities by adjusting the parameter 
^ of the RMF model that controls the quartic omega-meson 
self interactions ( Muller fc Serot|1996 1 and subsequently the 
high-density component of the EoS of SNM. It is set so that 
the SNM EoS is maximally stiff, resulting in maximum mass 
neutron stars of > 2.5Mq. The crustal EoS and crust-core 
transition densities are derived from the compressible liq¬ 
uid drop model (CLDM) for the crust ( [Newton et al.|[2013[ ) 
using the same nuclear matter model as the core EoS. EoS 
quantities obtained include the proton fraction Xp and mi¬ 
croscopic effective proton mass m* as a function of density, 
required to calculate the mutual friction parameter B from 
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equation |14| For a 1.4 Mq star, the softest EoS L =30 MeV 
produces effective proton masses that vary from 0.7 down 
to 0.3 from crust to core and the stiffest EoS L —120 MeV 
produces effective proton masses that vary from 0.85 down 
to 0.15 from crust to core. 

In Fig. 1^ we show the bulk properties of neutron stars 
for three representative values of the slope of the symme¬ 
try energy L = 30,60 and 90 MeV for the two families 
of EOSs: SklUFSU (maximum masses around 2 M©) and 
lUFSU (maximum masses around 2.5-2.6 Mq). In the left 
panel, the thickness of the crust ARc is displayed as a func¬ 
tion of mass M, while on the right the total radius of the 
star R is displayed as a function of mass M. 


4 RESULTS 

Before we get into discussions of the results, we remind read¬ 
ers that when we refer to the portions of the core coupled 
and uncoupled to the crust, we are referring specifically to 
the superfluid core neutrons. The small proton component 
of the core is assumed to always be tightly coupled to the 
crust. 

In Fig. 4 we show the averaged mutual friction profile 
B{f) as a function of cylindrical radius f throughout the core 
of a 1.4 Mq star (a) and a 1.8 Mq star (b). In each case we 
show the mutual friction profiles obtained with the L=30, 
60, 90 and 120MeV members of the SklUFSU EOS family. 
The shaded band covers the possible range of the threshold 
strength of mutual friction inferred from the 2000 and 2004 
Vela glitches. All regions of the star with a mutual friction 
greater than this threshold are coupled to the crust at the 
time of glitch, while those regions below the threshold are 
uncoupled from the crust at the time of glitch. 

As the slope of the symmetry energy L increases (gets 
stiffer), the strength of the mutual friction B{f) gets larger 
throughout the core of the star. A larger L leads to the sym¬ 
metry energy as a function density of iJaym(p) that increases 
more rapidly, leading to a larger proton fraction at a given 
density, which contributes to an increasing B{r) with den¬ 
sity (since B is almost proportional to Xp). In addition, the 
microscopic effective mass m* generically decreases with in¬ 
creasing density (from 0.7-0.85 at the crust down to 0.15-0.3 
in the core for a 1.4Mq star). Since B depends on rup — rtip 
and inversely on m*, this also contributes to an increas¬ 
ing B{f) with density. For the softest EOS shown, L — 30 
MeV, the symmetry energy begins decreasing with density 
just above saturation density, and therefore so does the pro¬ 
ton fraction. This behavior dominates the behavior of B at 
high density, which explains its non-monotonic behavior for 
L = 30 MeV. 

At the softer end of the range of L we consider, the 
mutual friction profile changes rapidly with variations in L. 
For the softest EOS shown, L =30 MeV, the non monotonic 
behaviour arises because the proton fraction decreases with 
increasing density from the crust to the core, while the effec¬ 
tive proton mass increases. For stiffer EOSs, both quantities 
increase from crust to core. 

For L = 30 MeV the predicted mutual friction pro¬ 
file lies almost completely within the inferred range. The 
upper end of the range of B{r) inferred from observation, 
1.3 X 10“^, implies that the core neutrons are entirely de- 


Table 2. Extension of the strong pinning region into the outer 
core required to satisfy the Vela constraint on G for a 1.4Mq 
neutron star for L = 30 MeV and different values of Yg. The 
extension is given as the baryon number density above the crust 
core transition density, n+. 


L(MeV) 

n+ {B = 1.3x10-4) 

n+ {B = 3x10“®) 

30 

0.016 

0.083 

60 

0.048 

0.051 

90 

0.032 

0.035 

120 

0.023 

0.025 


coupled from the crust at the time of glitch, and that the 
superfluid neutrons in the strong pinning region in the crust 
have only the rest of the crust and the core protons to spin 
up during the glitch event. At the lower end of the range, 
for a canonical 1.4 Mq neutron star a small portion of the 
outer core becomes coupled to the crust. For higher masses 
an increasing part of the inner core also becomes decoupled. 

If, on the other hand, we consider the lower end of the 
inferred range of B{f), 3 x 10“®, most of the core (within the 
inner ~ 10km) is coupled to the core. For L=60-120 MeV, 
most of the core has mutual friction significantly stronger 
than the inferred range, and hence a large fraction of core 
neutrons are coupled to the core at the time of the glitch, 
regardless of whether we consider the upper or lower bounds 
of the range. 

In Fig. 5 we plot the ratio of the moment of inertia of 
the coupled part of the core R (that which has a mutual 
friction parameter larger than the observationally inferred 
values) to the total moment of inertia of the star It as a 
function of mass M for the L=30,60,90 and 120 MeV mem¬ 
bers of the SklUFSU EOS. The shaded area in each plot 
spans the upper and lower bounds for R/It corresponding 
the fractions R obtained from considering the upper and 
lower bounds of the inferred range of B{r), 3 x 10“® as out¬ 
lined above in reference to Fig. 4. For Z/=30 MeV, the fact 
that the predicted mutual friction profile lies mostly within 
the inferred range of B{r) means that the ratio R is al¬ 
most completely unconstrained by the post-glitch evolution 
- as can be seen, between almost 0% and 100% of the core 
could be coupled to the crust during the glitch rise time. 
For intermediate-to-stiff symmetry energy behaviors at sat¬ 
uration density, L = 60 — 120 MeV, the range of R/R is 
inferred to be above 70% (above 80% for star > 1.4Mq). 

Given the moment of inertia of the predicted fraction of 
the core of the neutron star R coupled to the crust obtained 
as detailed above, we can calculate the ratio of the moment 
of inertia of the neutrons in the strong pinning region, as¬ 
sumed for now to be entirely within the inner crust, to the 
moment of inertia of the couple core fraction R/R. In Fig. 6 
we compare this value as a function of mass to the mini¬ 
mum required to explain the Vela glitch activity over the 
past 45 years G = 0.016, indicated by the horizontal dashed 
line in all four plots. It is important to note here that this 
comparison is meaningful only if the observed sizes of the 
glitch are a reasonable approximation to the actual change 
in frequency after the rise. This may not be the case if, as 
some simulations suggest ( Haskell et al.|20'l2 |, the rise time 
is significantly shorter than the current upper limit of « 1 
minute and a larger fraction of the core may be decoupled 
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Figure 5. Fractional moment of inertia of the component of the star coupled to the crust during the glitch Ic/hot as a function of stellar 
mass fori = 30MeV (a), L = 60MeV (b), L = 90MeV (c) and L = 120MeV (d). The ranges span the upper and lower bounds inferred 
from observation, and correspond to taking the upper and lower bounds of the shaded regions in figure]^ 


during a glitch. Furthermore the comparison is appropriate 
for the Vela pulsar, for which the activity is dominated by 
giant glitches where we can expect Mutual Friction to dom¬ 
inate the recovery. This is unlikely to be the true for other 
pulsars in which smaller glitches are frequent and are un¬ 
likely to lead to lag reversal. 

Returning to Fig. 6, the first thing to notice is that 
for no EOS in the SklUFSU family within the wide range 
of L considered does the predicted value of h/h exceed the 
required value to account for Vela glitches. Within this fam¬ 
ily of EOSs, the crustal neutrons are insufficient to explain 
the Vela glitches, using a core fraction consistent with the 
initially observed glitch recovery. 

Generically, the ratio h/h decreases with increasing 
mass as the crust becomes thinner while the radius stays 
relatively constant (refer to Fig. |^. At lower L the pre¬ 
dicted range for h/h spans an order of magnitude owing 
to the wide range of possible coupled core fractions. The 
upper bound of the range corresponds to the lower bound 
of h/h, which in the case of L = SOMeV is very small as 
the core is completely uncoupled over the glitch rise time. 
Because of this, the upper bound comes closest to matching 
the required level from the Vela pulsar glitches. In the three 
other cases, between 70 and 100% of the core is coupled, 
and the range predicted for h/h increases with L as the 


star’s crust increases; it does not rise enough to match the 
required level, though. 

If one allows the region where the neutrons experience 
strong pinning to penetrate the core to an arbitrary amount, 
la can be increased until the ratio h/Ic > 0.016. In Table 2 
we give the density above the crust-core transition density 
to which we must extend the inner boundary of the strong 
pinning region in order to satisfy the Vela constraint for 
a 1.4Mq star. Pinning could occur in the outer core, for 
example, on magnetic fiuxtubes. 

We next explore the extent to which changing the high 
density EOS changes the outcomes of our predictions. We 
now use the lUFSU family of EOSs with the parameter 
^ = 0.00 which leads to the stiffest EOS at super-saturation 
densities and a corresponding maximum mass of 2.6M©. The 
EOS and symmetry energy is the same as the SklUFSU up to 
saturation density. We show representative results in Fig. 7. 
In the left-hand panels we show the predicted ranges for the 
moment of inertia of the coupled core fractions to the total 
moment of inertia h/h, while in the right-hand panels we 
show the ratio of the moment of inertia of the neutrons in 
the strong pinning region to the moment of inertia of the 
couple core fraction h/Io, compared to the minimum value 
inferred from the Vela pulsar. We show results for a rep¬ 
resentative soft EOS L = 30MeV and representative stiff 
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Figure 6. Moment of inertia of the neutrons in the strong pinning region of the crust Is relative to that of the component of the 
star coupled to the crust during the glitch 7c as a function of stellar mass for L = 30MeV (a), L = 60MeV (b), L = 90MeV (c) and 
L = 120MeV (d). The ranges correspond to taking the ranges of 7c shown in figure The lower bound inferred from Vela’s glitch 
activity is shown by the horizontal dashed line. 


EOS L = 90MeV from the lUFSU ( = 0.00 family. The in¬ 
creased stiffness of the high-density EOS prevents the pro¬ 
ton fraction dropping and weakening the mutual friction for 
L — 30MeV as happens with the corresponding EOS from 
the SklUFSU family. The behavior of the coupled core frac¬ 
tion is thus similar to that of the stiffer EOSs including the 
L = 90MeV shown, above «80% for all masses. 


5 DISCUSSION AND CONCLUSIONS 

Evidence exists for a short timescale (~ 1 min) component in 
the post-glitch recovery of the 2000 and 2004 Vela glitches, 
suggestive of the recoupling of that portion of the stellar 
interior that was decoupled from the crust during the glitch. 
Under the assumption that the recoupling is mediated by the 
mutual friction force which results from electron scattering 
off the magnetized cores of superfluid vortices, we use the 
observed timescales and magnitudes to infer the moment of 
inertia of the region of the core coupled to the crust during 
the glitch. 

Within the framework of hydrodynamic models of vor¬ 
tex dynamics, which suggest the pinning region may be con¬ 
fined to the region of the inner crust where the vortices are 
completely immersed in the crustal lattice, we calculate the 


ratio of the moment of inertia of the pinning region of the 
crust to the moment of inertia of the core coupled to the 
crust during the glitch which we compare with the observed 
Vela glitch activity. We derive the neutron star EOS from 
the crust to the core, including the composition, microscopic 
effective proton mass and crust-core transition density, con¬ 
sistently using the same underlying nuclear matter EOS. We 
choose two families of EOSs: one derived from the SklUFSU 
parameterization of the phenomenological non-relativistic 
Skyrme model of the in-medium neutron-proton interactions 
and one derived from the lUFSU parameterization of the 
phenomenological relativistic mean-field (RMF) model. The 
Skyrme family of EOSs give maximum mass neutron stars 
of « 2.0Mq, while the RMF family give maximum mass 
neutron stars of « 2.6Mq. Each family is generated by sys¬ 
tematically varying the slope of the symmetry energy at sat¬ 
uration density over a conservative range L = 30— 120MeV, 
thus exploring the impact on our results of the current un¬ 
certainty in the EOS of neutron rich matter at saturation 
density. 

The main conclusion of our analysis is that taking the 
strong pinning region to be confined entirely to the crust, no 
EOS results in large enough glitches to explain the Vela pul¬ 
sar’s glitch activity for neutron stars above IMq. However, 
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Figure 7. Fractional moment of inertia of the component of the star coupled to the crust during the glitch 1^/hot (a,b) and moment 
of inertia of the neutrons in the strong pinning region of the crust h relative to that of the component of the star coupled to the crust 
during the glitch h (c,d) for the L = 30MeV (a,c), L = 90MeV (b,d) members of the lUFSU family of EOSs. The bands have the same 
meaning as figures 5 and 6. 


for the stiffest EOSs considered, neutron stars with masses 
slightly under IMq would marginally satisfy Vela’s glitch ac¬ 
tivity, and uncertainties on glitch sizes and rise times could 
lead to lower fractions of the core being coupled during a 
glitch and make low masses consistent with the data. 

We find that increasing the stiffness of the high density 
EOS, has only a small effect on the results for 1.4M0 neu¬ 
tron stars, generally increasing the moment of inertia cou¬ 
pled to the crust during the glitch slightly. This high density 
stiffening has two effects: increasing the radius of the star, 
and thus the moment of inertia of that part of the core that 
must be spun-up, and also increasing the crust thickness, 
and therefore the angular momentum reservoir which spins 
the star up. These two effect oppose each other, resulting in 
only a small effect overall. 

It would thus appear that even if one allows for part 
of the core superfluid to decouple during a glitch it is quite 
generally the case that, if standard mutual friction mecha¬ 
nisms are at work, strong entrainment in the crust limits the 
amount of angular momentum available for a glitch to below 
the level required to explain the glitch activity of the Vela 
pulsar. Allowing the strong pinning region to extend by up 
to 0.05fm“® into the core, would allow all EOSs to produce 
glitches consistent with Vela’s glitch activity. Physically this 


would correspond to the protons in the outer core being in 
a type-II superconducting state and the superfiuid neutron 


vortices being pinned to superconducting flux tubes (Link 
|2003[ [Haskell et al.||2013| , and contributing to the angular 
momentum reservoir available for a glitch. Note, however, 
that in this case vortex-flux tubes interactions could be the 
leading source of Mutual Friction in some areas of the core 
(Sidery & Alpar 2009 Haskell et al. 20141, thus invalidat¬ 
ing our assumption that Mutual Friction is due to electron 
scattering off vortex cores. 


The strength of mutual friction throughout the core in¬ 
creases with increasingly large symmetry energy slopes L. 
At low values of L (soft EOSs), the fraction of the core 
coupled to the crust at the time of the glitch is quite uncon¬ 
strained by the observations. At intermediate-high values of 
L (stiffer EOSs), the fraction of core coupled to the crust is 
constrained by observations to be above ~70% for realistic 
neutron star masses. 


The values extracted for the coupled fractions are con¬ 
tingent not only on our model for the crust-core coupling 
mechanism, but also the robustness of the observational 
data. As such, our results are tentative, as they depend on 
only two observations, one of which (the 2004) is particularly 
marginal. The fitting of exponential recovery terms can also 
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be called into question, as on short timescales the frequency 
evolution may be described by a very different functional 
form (Haskell et al. 2012 Link 20141. It is important to 
note that our conclusions also rely heavily on reliability of 
estimates of the frequency after a glitch. This may not be the 
case if the glitch rise time is signihcantly shorter than the 
current upper limit of « 1 minute, e.g. closer to timescales 
of a few seconds as simulations suggest ( Haskell et al.|20i2 1, 
and a significantly larger fraction of the core may be de¬ 
coupled during a glitch. If we allow all the core neutrons 
to be entirely decoupled during the glitch (requiring cou¬ 
pling timescales of less than 5-lOs depending on the EoS), 
we still do not obtain consistency with a IAMq Vela pulsar 
except for the stiffest EOS, L =120 MeV. However, given a 
reasonable uncertainty of order 0.01 fm“^ in the crust-core 
transition density (Newton et al. 20131 at a given value of L, 
our model would be consistent with Vela being a M = 1.40 
neutron star if L > 70 MeV and the coupling timescale is 
sufficiently fast to decouple > 70% of the core neutrons. 
From our calculations, this requires the coupling timescale 
to be « 10s or less, as opposed to the minute timescales 
estimated from the observed glitch recoveries. Consistency 
with Vela being a M = 1.0© neutron star can be achieved 
for L > 50 MeV, provided the coupling timescale is « 30s 
or less, enough to decouple > 80% of the core neutrons. 

The one aspect of our analysis that is not consistent 
with the underlying EOS is the use of equation for the 
mutual friction coefficient; although the proton fraction and 
proton effective mass are derived from the underlying EOS, 
the form of the equation is a ht to a single microscopic calcu¬ 
lation. Although we expect the basic functional dependence 
on quantities not to change much, the detailed form may de¬ 
pend on the EOS. One should also note that the values for 
the mesoscopic effective neutron masses in the crust, which 
encode the entrainment effect, are also not consistent with 
the EOS used, and could substantially alter the estimates 
of the angular momentum available in the crustal super¬ 
fluid. Another important assumption of our model is that 
the protons and electrons in the core are coupled to the crust 
on short timescales by the magnetic field. This assumption, 
however, depends strongly on the magnetic field geometry 


and may not hold for all regions of the core (Glampedakis 
|fc Lask^|2015[ ). 

Nevertheless, it is likely that the post-glitch response 
of glitches bears the imprint of the crust-core coupling dy¬ 
namics, and our study illustrates one way in which such 
information can be extracted. 
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